Momentum-dependence of charmonium spectral functions from lattice QCD 
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We compute correlators and spectral functions for ,J/tp and 77c mesons at nonzero momentum 
on anisotropic lattices with Nf = 2. We find no evidence of significant momentum dependence 
at the current level of precision. In the pseudoscalar channel, the ground state appears to survive 
up to T ~ 450MeV or 2.1T'c- In the vector channel, medium modifications may occur at lower 
temperatures. 



I. INTRODUCTION 

The survival and properties of charmonium bound 
states in the quark-gluon plasma is a hot topic in the 
interpretation of experimental results from relativistic 
heavy-ion collisions. A variety of models have been put 
forward to explain the observed pattern of suppression 
of the J/ip yield, including sequential suppression [l|, 
charmonium regeneration through coalescence of uncor- 
related cc pairs 0t1|; and combinations of these two 
mechanisms @. 

Results from lattice simulations ^-Q suggest (though 
far from decisively [l^l) that S-wave ground states (J/V' 
and rjc) survive up to temperatures close to twice the 
pscudocritical temperature T^, while excited states {tp') 
and P-wavcs (xc) niclt at temperatures close to Tc, in 
rough agreement with the sequential suppression sce- 
nario. These results are however not sufhcient to rule 
in favour of one or other scenario. 

An additional handle on the problem may be gained 
by considering charmonium states which are moving with 
respect to the thermal medium. The two main scenarios 
(regeneration and sequential suppression) predict differ- 
ent transverse momentum and rapidity dependence of the 
charmonium yield. This is primarily due to the momen- 
tum dependence of the reaction rates included in these 
scenarios, but in order to complete the story it will also 
be necessary to investigate any momentum dependence 
of the baseline survival probability of charmonium bound 
states. This paper provides the first steps in this direc- 
tion. 

The analytical structure of QCD correlators implies 
that both spatial and temporal correlators are governed 
by the same spectral function. However, in order to de- 
scribe spatial correlators and screening masses, it is nec- 
essary to also include the spatial-momentum dependence 
of the spectral function. 

Finally, computing spectral functions at nonzero mo- 
mentum may also facilitate the determination of trans- 
port coefficients such as the heavy quark diffusion con- 
stant |ll| , as the zero- frequency transport peak is broad- 
ened and shifts to nonzero frequency. 

We will work on a spatially isotropic lattice, which 
means that the anisotropy of the plasma created in 
heavy-ion collisions is not taken into account. Our mo- 



mentum is best interpreted as the transverse momentum, 
and we will therefore not be able to say anything at this 
point about the rapidity dependence of any obscrvables. 



II. FORMULATION 

We have simulated two degenerate flavours of 0{a)- 
improved Wilson fermions, with ra-^jrap = 0.54 (corre- 
sponding approximately to the strange quark mass) , and 
with spatial lattice spacing a^ = 0.162fm and anisotropy 
^ ~ asja-r = 6. The lattice action is described in more 
detail in Ref. [12|, and the high-temperature simulation 
parameters are discussed in Ref. 11311 . Our parameters 
here correspond to Run 7 in Ref. [13| . The same fermion 
action is used for both light sea and heavy valence quarks. 
The lattices, corresponding temperatures and number of 
configurations A'cfg used in this study a re g iven in TablelD 
The main difference compared to Ref [l3| is that all our 
data arc now obtain using the tuned (Run 7) parameters 
and on a spatial volume of 12"^ in lattice units. 
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TABLE I: Lattice volumes, temperatures and number of con- 
figurations used in this simulation. The separation between 
configurations is 10 HMC trajectories, except for the Nt — 80 
runs where configurations were separated by 5 trajectories. 
The lattice spacings are as ~ 0.162fm, Ut — 7.35 ± 0.03GeV. 

Spectral functions were computed using the maxi- 
mum entropy method as described in |13| . with the 
modified kernel introduced in Ref. [ij]. The charmo- 
nium correlators were computed using all-to-all propaga- 
tors with dilution in time, color and varying levels in 
space and spin. The momentum values are given by 
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FIG. 1: The mass splittings of the low- lying charmonium 
spectrum at zero temperature, with the known experimental 
values and the DD threshold also indicated. 
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TABLE II: Operator basis used to obtain the T — charmo- 
nium spectrum. The definitions of Si,pi and ti are given in 
Ref. \m. 



p^ = f ^'^ \ n^, with ri^ = 0,1,2,3,4, corresponding 

to p = 0, 0.66, 0.93, 1.14, 1.32 GeV. We consider only the 
pseudoscalar (77c) and vector ( J/V") channels. Our corre- 
lators have not been renormalised, so the vertical scale 
on the figures presented is arbitrary, but the shape of the 
curves is not. 



III. RESULTS 
A. Zero-temperature spectrum 

In order to establish a baseline for our high- 
temperature results, we first present results for the spec- 
trum and dispersion relation at zero temperature. Taking 
advantage of the all-to-all propagators, we used the varia- 
tional method [1^, [la| with a large basis of interpolating 
operators to compute the low- lying spectrum of S-, P- 
and D- waves. This is shown in Fig. [TJ along with the ex- 
perimentally known masses. The operator basis used in 
the spectrum study is summarized in Table |lTl Taking 
the difference between the spin-averaged IS masses and 
the ^Pi (/ic) mass to set the scale, we find the inverse lat- 
tice temporal lattice spacing to be a~^ = 7.35f3)GeV|19l|. 
In Figure [l]the highest lying radial excitations in each 
channel are colored red to indicate that they are contam- 
inated by higher excited states. The difference between 
our results and experiment is largely due to the fairly 
coarse spatial lattice spacing used. Including a clover 
term in the action has been found to improve the agree- 
ment with experiment jlTJ . 

Figure [2] shows the ground-state pseudoscalar (77c) dis- 
persion relation. The bare anisotropy for the valence 
charm quarks is ^g = 5.90, which on the 8'^ x 80 lattice 
had given a renormalised anisotropy of 6. However, on 
the larger (12^ x 80) lattice we find that the dispersion re- 
lation gives us a renormalised anisotropy of 5.829±0.05, 
giving rise to a 3% systematic error in our final results. 
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FIG. 2: The r\c dispersion relation. 



B. High-temperature spectral functions 

In order to facilitate comparisons between spectral 
functions at different temperatures and momenta, we 
have determined the spectral functions p{ui,p;T) at all 
T and p using the same default model, m{uj) = mow^ 
with ariTiQ = 6 and aT-rriQ = 4 for the pseudoscalar and 
vector channel respectively. These values are close to 
the 'best' values determined from a one-parameter fit of 
the correlator for each temperature and momentum con- 
sidered. We have repeated the analysis for a range of 
different default models to study the robustness of our 
results. 

As in Ref. [ij], the first two timeslices were discarded 
in the analysis since these will be dominated by lattice 
artefacts. The effect of varying the time range has also 
been investigated. 

Figure |3] shows the spectral function in both channels 
at zero momentum for the various temperatures. These 
results may be compared to Figs. 5 and 6 in Ref. [13|. 
In the pseudoscalar channel, our results suggest that the 
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FIG. 3: The pseudoscalar (top) and vector(bottoin) spectral 
function at zero momentum for various temperatures. 



ground state {rjc) survives up to the highest tempera- 
tures accessible in our simulations, while the J /ip (vec- 
tor state) appears to dissolve for T > 370MeV or l.lTc- 
It must be noted, hovifever, that the uncertainty in the 
MEM procedure increases as the temperature increases, 
and it is therefore not possible at present to say for cer- 
tain whether the pattern observed at the highest tem- 
peratures is genuine or merely reflects the inability of 
the MEM algorithm to determine the spectral function 
given the available data. This remains the case for all 
momenta and in all channels. 

Figure 0] shows the pseudoscalar spectral function for 
different momenta at T = 230 MeV. Using Tmin = 3 
(lower panel) , we see a clear structure of peaks ordered 
by momentum, with the separation between the peaks 
corresponding reasonably well to the zero-temperature 
energy levels. The peak positions, however, appear to 
be shifted compared to the zero-temperature energies. 
This is most likely due to the maximum entropy method 
not being able to resolve the full detail of the spectral 
function for the available data. This is supported by the 
upper panel of Fig. HI showing the spectral functions de- 
termined using Tmin = 2. Here the peak position for the 
lowest non-zero momentum corresponds precisely to the 
zero-temperature energy, while for the other momenta 
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FIG. 4: The pseudoscalar spectral function at nonzero mo- 
mentum for T = 230 MeV, using timeslices 2-16 (top) and 
3-16 (bottom). The zero-temperature energy levels are also 
shown for comparison. 



the peaks are smeared out and shifted to higher energies. 

Figure [5] shows the pseudoscalar spectral function for 
different momenta at two higher temperatures. There is 
still some evidence of a surviving ground state peak at all 
momenta, but uncertainties in the MEM reconstruction 
means that we are no longer able to resolve the ordering 
of the states given our current precision. 

At nonzero momentum, the vector meson correlator 
is decomposed into transverse and longitudinal polarisa- 
tions. 
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The transversely and longitudinally polarised J/V' may 
in principle behave differently in the medium. We have 
therefore analysed the two separately. Figures [HI and 
[7] show the transverse and longitudinal spectral func- 
tions for our largest and smallest momentum, respec- 
tively. The longitudinal spectral function for the lowest 
momentum {p = 0.66GeV) shows clear evidence of a sur- 
viving peak corresponding to the J /ip ground state at 
the two lowest temperatures, but there are indications 
of medium modifications already at T = SOOMeV. The 
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FIG. 5: The pseudoscalar spectral function at nonzero mo- 
mentum for T = 368 MeV (top) and T = 459 MeV (bottom). 
The zero-temperature energy levels are also shown for com- 
parison. 



transverse spectral function, on the other hand, appears 
to be subject to medium modifications (indicated by the 
softer ground state peak) already at the lowest tempera- 
ture, just above T^. 

It is viforth emphasising that the transverse spectral 
function at T = 230MeV, p = 0.66GeV displays little 
or no sensitivity to the choice of model function, so we 
consider our result in this case to be robust. At higher 
momenta and temperatures, the MEM reconstruction 
of spectral functions is no longer stable with regard to 
variations in the model function and time range used. 
Our results for both transverse and longitudinal J/ip at 
p = 1.32GeV should therefore be taken as no more than 
indicative, even at the lowest temperature. With this 
proviso, the indications are that the longitudinal J/ip 
survives at least up to 265MeV (1.2Tc) and is most likely 
melted by T = 400MeV, while the transversely polarised 



J/ip may experience medium modifications already close 
toT,. 

IV. DISCUSSION 

We have presented first results for charmonium spec- 
tral functions at nonzero momentum. Our results suggest 
that the IS pseudoscalar meson (r]c) survives up to tem- 
peratures close to twice the pseudocritical temperature 
of QCD, for all momenta. No substantial momentum 
dependence was found. 

In the vector channel, there appears to be a distinction 
between the transverse and longitudinal channels, with 
the longitudinally polarised J/ip experiencing smaller 
medium modifications. Again, no substantial momen- 
tum dependence was found, although the reconstruction 
of the spectral function became progressively more un- 
certain with increasing momentum. Therefore, some ad- 
ditional momentum dependence can not be ruled out. 

Whether the difference between transverse and longi- 
tudinal polarisations is real or merely a reflection of the 
uncertainty in the MEM procedure given the data used 
in this study, still needs to be determined. 

We are in the process of generating configurations with 
smaller lattice spacing. This will provide greater tem- 
poral resolution, leading to a more reliable determina- 
tion of spectral functions from imaginary-time correla- 
tors. The finer lattice is also expected to bring the zero- 
temperature spectrum in closer agreement with exper- 
iment, and allow a clearer separation between physical 
features and lattice artefacts in our spectral functions. It 
will also allow us to access higher temperatures, where 
all model studies up to now have found the charmonium 
ground states to be dissolved. 

We are also computing correlators of the conserved vec- 
tor current, which will provide the correct nonpcrturba- 
tive renormalisation of the vector operator used in this 
study, and permit a quantitative determination of the 
charm quark diffusion rate and the charmonium contri- 
bution to the dilepton rate. 
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FIG. 6: The transverse (left) and longitudinal (right) vector meson spectral function at p = 0.66GeV for various temperatures. 
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FIG. 7: The transverse (left) and longitudinal (right) vector meson spectral function at p = 1.32GeV for various temperatures. 
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